Diagrammatic method for investigating universal behavior of impurity systems 
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The universal behavior of magnetic impurities in a metal is proved with the help of skeleton 
diagrams. The energy scales are derived from the structure of the skeleton diagrams. A minimal set 
of skeleton diagrams is sorted out that scales exactly. For example, the non-crossing approximation 
for the Anderson impurity model can describe the crossover phenomenon. The universal Wilson- 
number is calculated within the non-crossing approximation. The method allows for an assessment 
of various approximations for impurity Hamiltonians. 



I. INTRODUCTION 



Magnetic impurities in metals show universal behavior at low energiesu. (The Hamiltonian of such systems consists 
of at least one conduction band of width D to which the impurity couplesa. At temperatures ksT <C D, the scaled 
observable is independent of details of the host system such as its band structure. 

A long standing question in the physics of such systems is: How can the observed universal behavior directly be 
established from the original model Hamiltonian such as the Anderson-impurity modelcl, together with a reasonable 
accurate description of observables? Hitherto the original Hamiltonian has been replaced by another one which is 
' O i more accessible: 

(a) Within the Bethe ansatz the impurity part of the system's thermodynamics can be derived. However, the 
original model has to be replaced by one with linear dispersion. The spectrum of eigenvalues is cut off at D' which 
is in general not identical with the band width D of the metallic host in the original modeHi A relation between D 
and D' has so far been established only for the Anderson impurity model! 

(b) Within the numerical renormalization groupfl, the impurity is coupled to a half-infinite chain with hopping 
matrix elements vanishing as A _n , n = 1, 2, . . . , A > 1. In the limit of A = 1 one would recover the original model. 

■ However, the numerical results have to be extrapolated to that limit because the length of the chain to be diagonalized 
' numerically would eventually become too largeaH. 
t~^. . Hence both methods do not prove that the observables of the original model behave universally. On the other hand, 
there are the diagrammatic approaches to the impurity problem!! They allow us to construct approximations for all 
observables of the original model. The Dyson equation 



R~\z) =i? fl - 1 (z)-S(z) 



is invariant under a certain rescaling of the propagators R and Rq, self-energies S, and coupling constants. With the 
help of the diagrammatic renormalization group, perturbative results for the propagators, self-energies, and vertex 
parts are then fitted to the Dyson equation. In that way, scaling laws for the propagators result, from which the 
universal behavior in the perturbative high-energy regime of the modelcl follows. This method amounts to summing a 
certain subclass of diagrams, with naked propagators. However, the procedure breaks down when perturbation theory 
fails. 

For the nonperturbative region dressed propagators are necessary. This requires the use of skeleton diagrams, oth- 
erwise the definition of a self-energy itself would become ambiguous. However, the scaling of the dressed propagators 
^ | is unknown, precisely because beforehand the Dyson equation would have to be solved. _ _ 

In this paper, this difficulty is overcome by utilizing the variational principle of Luttinger and Warcu, KuramotoE3, 
and Baymtll, by which any observable can be expressed in terms of skeleton diagrams. This is exemplified at the 
Anderson-impurity Hamiltonian. 

It turns out that the skeleton diagrams of the second order are already sufficient to calculate the exact energy 
scales. This proves universal behavior for this model. The second-order skeleton diagrams are therefore the minimal 
class of diagrams which have to be summed, in order to describe the crossover. 
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II. HAMILTONIAN AND DIAGRAM TECHNIQUE 



As the standard model describing magnetic impurities in metals, the Anderson-impurity HamiltonianEl is considered, 
with a half-filled conduction band of constant density of states p and infinite Coulomb repulsion at the impurity site: 

H = H C + H f + H 1 , 

He — ^ t j e p c pm c pm 
|e p |<£>,m 

H f = e f^2 fmfm 
m 

Hi = ^J2>(4 m b + f m + H.c). (1) 

c pm creates a conduction electron with internal quantum number m = 1 . . . N, momentum p, and energy e p which is 
cut off at ±D. \m) — /+|vac) denotes a magnetic configuration of the impurity and |0) = 6 + Uac) the non-magnetic 
one, their energy-difference being ef. is a fermionic operator and b + the Coleman bosor£3. F£ = f+b creates 
an electron at the impurity site. Double occupancy at the impurity site is suppressed by imposing the constraint 
715 + rif = 1. The impurity hybridization with the conduction band is proportional to V. The Boltzmann constant is 
set to unity so that temperature is measured in units of energy. 

In order to simplify the subsequent derivations, the conduction band of model (|l|) is assiHaed to have a constant 
density of states p with a sharp cutoff at energie s ±Z ), as the other approaches to scaling dcBEffi. Universal behavior 



should not depend on this assumption!. In Sec. VII it will be shown that the energy scales of the system indeed do 
not change as long as the density of states is finite at the Fermi energy, and is sufficient structureless to have only one 
energy scale D. 

The principal object of concern is the resolvent Rf(z) = (Q(z)) c where Q is defined as 

Q(z) = \/{z-L c -H f -Hi). 

Here the superoperator L c acts on an operator X of the Hilbert space as L C X = [H C ,X], and () c indicates the 
thermodynamic average with respect to H c . The propagators for the occupied and unoccupied ionic configurations 
are 



With the help of the identity^ 



R m (z) := (m\Rf(z)\m) 
R (z) := (0\R f (z)\0). 



e -0H _ e -P{L a +Ht+Hx) e -PH e (2) 



the impurity part Zf of the partition function can then be represented as a line integral, the path of integration 
encircling all poles of the integrand: 

-■- TT ^ m -Tr J ^R f (z). (3) 



1 ' Tr c e-^ " 1 J 2ni 

The well-known diagrammatic techniqueEl follows if Q is expanded in a geometric series in V, then L c acting on the 
conduction-electron operators in Hi is evaluated to give the energy denominators 



and finally Wick's theorem applied to evaluate the thermodynamic average with respect to H c . This can be casted 
in a diagrammatic language. The naked propagators and vertices are shown in Fig. |l| Because the impurity site can 
alternatively be empty or singly occupied, every diagram has a spine of alternating b and / propagators. Consequently, 
all diagrams where a conduction-electron propagator would have a self-energy are excludedEfl, because they do not 
fulfill t he con straint n^ + nf — 1. Alternative approaches where this constraint is not exactly enforced will be discussed 



in Sec. VIII 



Within the variational principleuEil, a functional T of the dressed one-particle propagators is defined. Because 
the conduction-electron propagators carry no self-energy, the variational principle has been adapted to the present 
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diagrammatic technique^ such that T becomes a functional of the dressed one-particle impurity-propagators only. 
At the saddle point with resppet to variations of the i?o,mj the functional T equals Zf, and the Dyson equation holds 
as a self-consistency equationl!3 



R f {z)- 1 =z-Hf-V f {z). 
For the impurity-part of the partition function, the functional is given by 

dz 



(4) 



T = (3 Tr 



2m 



"El 1 -; 



W(z)R f (z) 



In 



-ff/-£4 n) W 



(•5) 



Here denotes all nth order self-energy diagrams of Rf, expressed in terms of skeleton diagrams. Skeleton diagrams 
are diagrams where all self-energy insertions have been removed. Hence a skeleton diagram becomes a functional of 
the dressed propagators. 

The variational principle can be interpreted in the following fashion: If T depends on some parameter A suchj-as 
the hybridization then, at the saddle point, it depends on A only explicitly and not implicitly via the propagatorsE3: 



~dX 



<9T 
dA 



ST dR , 
5Ra, m d\ 



<9T 
~d\' 



Explicitly, from Eq. (0) it follows at the saddle point that 



dZj_ 
dX 



dX 



where the A dependence of Rf can be discarded. 

In order to study the dependence of Zf on a parameter A, it therefore suffices to consider all skeleton diagrams of 

the type S^-Rcm- A skeleton diagram for £Q 2 ^i? , i.e. of second order, is shown in Fig. |^. 

Approximations fulfilling the variational principle can be generated by using a subclass, a so-called family of skeleton 
diagramsEl, which contains with each skeleton diagram for EQ^-Ro.m all others with cyclic permuted vertices as well. 

For example, all skeleton diagrams of a given order Sq" '^Rq m form such a family. 

For instance, if only skeleton diagrams of second order are kept in Eq. ([5]), this amounts to summing all diagrams 
with bare propagators and noncrossing conduction-electron lines and is called the NCAE3. The self-energies of the 
NCA, E$(z) = S^ ) [i? (z)] and s[, 2) (z) = E|, 2) [i?i(z), . . . , R n (z)}, are then (/ denotes the Fermi function) 



V 2 P 
N 

V 2 P 
N 



D 



f(e)Ro(z + e)de, 



(6) 



-D 
N 



V / f(e)R m (z + e)de. 

„,_i J-D 



m — 1 ' 



The functional T^ 2 ^ has then the form 

T (2) 



In 



*-£< a) (z)]+5> [z-e f 



To construct a variational principle for other observables the Hamiltonian has to be coupled to external fields suitably 
chosenEI This will be discussed in Sec. |v|. 
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III. ENERGY SCALES 



A. First scaling equation 

T as defined in Eq. ^ depends explicitly on e/ i^ply via the term H j = £/ Yl m fmf m - Consequently, one has for 
the impurity part Ff = —TlnZf of the free energyli-3: 

Z,e f d ef F, = Tr f H f j> ^e^Rf(z). (7) 
T depends explicitly on V vi the prefactor V 2n of the 2nth order self-energy : 

Z fP V% v iF f = ~Tr f j ^-e-^ f {z)R f {z). (8) 

To determine the T dependence, the internal integration variables z and e as in Eqs. (|^), and (|6|) are replaced by Tz 
and Te, respectively. The variational principle remains unaffected. T depends now explicitly on T via the prefactor 
T n of a skeleton diagram of 2nth order, the term Tz in the logarithm, and the integration boundaries as in Eq. (^) 
change to ±D/T. Therefore 

Z f (F f -T±F,) =^,/^e-^E /W fl /W 
8 f r\7 

+ZfD dD Ff + TTfHf f 2ri e ~ PZRf{z) - (9) 

Inserting Eq. (0) and Eq. |j) into Eq. @ yields 

F '-( T w + ^^ +e 'i- f +D w) F '- (10) 

This equation expresses the scaling of Ff with respect to the energy scales D, e/, and pV 2 : 

F(T,p,V,e f ,D)=Tf' 
This is henceforth called the first scaling equation. 



pV 2 ' e f ' D 



B. Second scaling equation 



The central issue of this work is to prove and describe the universal behavior of impurity systems. To show the 
universal behavior for Ff, the functional T is examined for largej-but finite D, that is, when D becomes larger than 
all other energy scales of the system, the so-called universal limitd. 

T depends explicitly on the cutoff D only via the integration boundaries ±D of the integration over the conduction- 
electron energies, as can be seen in Eq. (p|) for the second-order skeleton diagrams. 

At first only those second-order self-energies are kept in T. With the help of the spectral densities po, m of i?o,m 
one has 

D dD F i 2) = wi^ I e ^ duJ D fO>D)M») 

■m—1 

x $R (lu + pD) + po(u) ^R m {uo + pD)] (11) 

where 5Ri?o,m denotes the real part of i?o,m- I n the universal limit, in particular T/D — > 0, so that only the terms oc 
/(— D) m 1 survive. Furthermore, in this limit the weighted spectral densities e _ ^ w po^ m (uj)/Zf contribute significantly 
only for frequencies less then the impurity part of the ground-state energy Eq < 0, because e^ E °Zf tends to 1 for low 
temperatures. From perturbation theory it follows^ that e~^" Pf(w) [Zf vanishes asymptotically as 1/uj 2 for large, 
negative u). Hence it contributes significantly to the integral in Eq. (O) only for 
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-VD <lo<E . 

In this frequency interval, 3lRf(u) — D) can be replaced by its bare counterpart l/(u> — Hf — D) as {—1)/D: 

D dD F f 2) = ~~NZ~~ ^ / e ^ dU [PnM + PoH] ■ 

f m=l J 

Together with Eq. (0), the following scaling relation is obtained: 

= PV 2 (1 - 1/N) ^-F? - P V\ (12) 

This is called henceforth the second scaling equation. Next it is shown that in the universal limit all families of 
higher order skeleton diagrams of T are irrelevant, by which is meant here that their contribution to the logarithmic 
derivative in Eq. (|l2] ) vanishes as 0(1/ D). For the proof see Appendix [a|, the families of skeleton diagrams being 
needed to enforce the variational principle. 

The result can be made plausible by casting it into the language of diagrams, as in Fig. [| Differentiating logarith- 
mically with respect to D means removing one curved conduction-electron line and replacing the internal propagator 
by its value at the cutoff cx 1/D. Therefore this diagram contributes cx D/D to the logarithmic derivative. 

For a diagram of higher order than two such as in Fig. |3j there lie under each conduction-electron line at least 
three propagators, because otherwise this diagram would have a self-energy insertion. Hence its contribution to the 
logarithmic derivative is cx D/D 3 and can be neglected for large D. . 

Equation (|l^) therefore holds as well for the exact Ff. It describes the Haldane scalinglJ for the difference Ff — Eq, 
that is, it depends on D and e/ only via 

(F f -E )(T,V,p,e f ,D) = (F f - E )(T, V, p, e}), 



* f = e f + (l- l/N) P V z \nD/{pV z ). 
By inserting this into Eq. (|l0|), in the magnetic limit — e/ 3> pV 2 the scaling law 

F f - E = Tg(T, V, e f , D, p) = Tg(^, I-) (13) 

follows, with the quantity T — irpV 2 /N and (up to a numerical factor) the Kondo temperature^ 



T K - D VpV 2 /D exp [e f /pV 2 ] . (14) 

This proves the universal behavior for the free energy of the Anderson Hamiltonian as well as that the energy scales 
r and Tk are the exact ones. In particular, it has been shown that the NCA preserves the exact energy scales of the 
system. Because the NCA is exact up to orders V 4 and 1/N, an inclusion of families of skeleton diagrams of higher 
order in Eq. ( |l3| ) will slightly change g but not the energy scalesn-and consequently will not alter the approximation 
qualitatively. Such an extension of the NCA has been performedtSl and within the errors of the numerical calculation 
the scaling law ([l3]) as well as a slight change in the respective function g have been verified. 

There holds an analogy to the usual diagrammatic renormalization group, as described in Sec. Q. It turns out that 
from a certain order in perturbation theory on, the energy scales obtained by that method do not change any more in 
the universal limit, while, for an observable the form of the respective function g can still change. However, the point 
is that this scheme can only be used for high temperatures T ^> Tk, below which perturbation theory breaks down. 
Hence a necessary condition for a diagrammatic technique to describe the universal behavior of impurity systems is 
that it includes families of skeleton diagrams including those of second order. This result can straightforwardly be 
extended to the case of finite magnetic fields. 

IV. SCALING WITH MAGNETIC FIELD 

The influence of a magnetic field h on the host metal is of order h/D cx 1/D, and can be neglected. Hence, to the 
Hamiltonian of Eq. ([!]) only 
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Hf — > H f + g[i B h^2 m f+f m (15) 

in 

is added, where g denotes the g factor of the / electron and pb is the Bohr magneton. The functional T will explicitly 
depend on h only via the new term in Hf. Consequently the first scaling equation ( p^| ) changes to 

F = T ^r F + py2 ^ F + ef lh F 

The second scaling equation ( p"2| ) remains unaffected. The scaling law for the magnetization M{h) = —(d/dh)Ff 
therefore reads 

M(T, V, e f , D, p, h) = M (J, . (17) 



V. PROOF OF UNIVERSALITY 



To obtain the scaling of any observable, the Hamiltonian H has to be coupled to suitable external fields. However, 
only the skeleton diagram of Fig. || is relevant for large D. Analogous scaling equations such as Eqs. ( p^| ) and ( |l3| ) 
can thus be derived for any other observable, hence proving universality for the Anderson model. 

As an example, take the / propagators as a function of imaginary time t, 

G fm (r) = -(T T F m (r)F+(0)). (18) 

Here T T is the time-ordering operator, and () denotes the thermodynamic average with respect to H. With the help 
of Q(z) one can reproduce the well-known integral. representation for G/ m which coincides for fermionic Matsubara 
frequencies with the Fourier transform of G / m (t)Ej, 

G fm (co) = fe-P'^{{0\Q(z)\0){m\Q(z + u)\m)) e . 

To obtain a variational functional for Gf m , H is coupled to an auxiliary fermionic field 

H + lu^+V + V^/pX(^ + b + f m + H.c). 

The / propagator is obtained from all diagrams of second order in A that is, after removing the ^ propagatorEI. A 
functional can be constructed analogously to Eq. (||), to give, at the saddle point, 

Tr/TrcTrvfre^ 

Gf follows as 

d 
8X* 



F*=pV 2 f{u)G fm (uj). (19) 



A=0 

The scaling equations for Fq,, and via Eq. ( |l9| ) for Gf, can be established in the same manner as for Ff. In particular, 
Eq. ( |i~(i| ) now reads 

T ^f + ^ + ^ + py2 a^ + D w + = °- (20) 

In the universal limit of large D again only the skeleton diagrams of second order are relevant, and the analog to 
Eq. (||) holds for G f , too: 

D^G fm =pV 2 (l-l/N)£-G fm . 

In the magnetic limit — e/ 3> T where T 3> Tk it follows that the Abrikosov-Suhl resonance and the impurity part 
of the resistivity computed within the NCA scale with the exact Kondo temperature. This was numerically observed 
in Ref. There, the energy scale was determined as the maximum value of the spectral density of G/ m (w) in the 
magnetic limit of the Anderson model, and assumed to be proportional to Tk- Here this has been proved. 
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VI. DESCRIPTION OF CROSSOVER: WILSON NUMBER 



Within this scaling method the crossover phenomenon!] can now be described entirely within a diagrammatic 
approach. Even for the skeleton diagrams of second order an analytical solution of the self-consistency equations 
does not seem possible for finite-|teinperature. However, for zero temperature the well-known expressions for the 
ground-state energy of the NCALL9EII for the Anderson model can be evaluated analytically (see Appendix |b|) in the 
magnetic limit. In this limit, it follows from Eq. ( |l7j ) that for zero temperature the scaling law 



M(V,e f ,D,p,h)=M 



h 

Tk 



(21) 



holds, which holds as well for the NCA, as shown above. For low magnetic fields, the law-field energy scale Tj, can 
be fixed unambiguously by the static magnetic susceptibility at vanishing magnetic fieldcl, 

X(0) = \$/T L . 



With the help of the analytical NCA results of Appendix ||, an analytical expression for Tl is obtained: 

T K /T(1 - 1/N), 



rpNCA 



(22) 



with r being the gamma function. The result coincides up to order 1/N with the result Tl = TrT(1 + 1/N) which 
is believed to be exactS. 

For high magnetic fields gpsh ^> Tk, the NCA reproduces the perturbation theory up to V 4 . With the definitions 
N = 2j + 1 and J = V 2 /\ef\, for the impurity part of the magnetization the expansion 



M(h) 



i-El 

N 



( P J) 2 ln g» B h 



N 



+ 



{pjf 

N 



N-l 



N(N 



1) ^ 



mln(m) 



0(J 3 ) 



is obtained. The exact scaling (|21j) of the NCA renders it possible to -reexpress in the perturbation expansion J in 
terms of h/Tx- The high-field energy scale T# is fixed unambiguouslyH by requiring that terms cx 1/ln 2 [gpBh/Tn] 
to be absent in the resulting asymptotic expansion. One arrives at the well-known asymptotic renormalization-group 
result for the magnetization in high magnetic fields 



M(h) 



In In 



= 1 - 



-O 



gp,Bh 
T H 



N , gpBh 
m ■ 



N 2 



T H 



In" 



gpBh 
T H 



lnln^\ 
T H 

, ln 3 ^ , 
V T H ) 



(23) 



and Th is given by 



T H = T K exp 



N-l 



-1/N - 



N(N-l)^ 



— to ln(m) 



(24) 



The Wilson ratio W = Th /Tl characterizes the crossover from the high- field region where the impurity reacts as an 
asymptotically free magnetic moment, to the low-energy region where the impurity is screened. Hence the Wilson 
ratio coincides with the exact result up to order 1/N: 



W 



NCA 



1 



ir 



1/N) 



r(i- i/jv)r(i- 

= 1 + (9(1/A 2 ). (25) 
By a Schrieffer- Wolff transformationQ, this result can be extended to the Coqblin-Schrieffer model (see Appendix [(]). 
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VII. HOW MUCH DOES SCALING DEPEND ON BAND STRUCTURE? 



The universal behavior of impurity models like the Anderson model should not, depend on details of the host's 
density of states as long as its band has a finite density of states at the Fermi energvu. This was realized for Jjupdel (|l|) 
by assuming a constant density of states with a sharp cutoff at ±£>, as the other approaches to scaling dc&lffl. There 
it is assumed that the energy scales of the system do not change as long as the density of states is finite at the Fermi 
energy, and scales as p(e) — p(e/D), which means that the band is sufficient structureless to have only one energy 
scale D. 

Here this is proved: The first scaling equation ([To]) can be taken over because of 

d d 
T—p(eT/D) = -D—p(eT/D). 



For the second scaling equation (12), the term 

Y:=V 2 I p(e/D)A 2 (u + e)f(e)de 



as in Eq. ( |A2| ) in Appendix ^ is examined, which represents an integration over a conduction-electron line. A2 stands 
for the real part of the m propagators which lie under this conduction-electron line. Differentiating logarithmically 
with respect to D and substituting e :— uD yields 

DgjjY = (-D) / up'{u)f(uD)A 2 (uj + uD) du. 



The integrand can be neglected outside the interval —\JD < w < as in Sec. [II B . If p is sufficiently smooth 
around the Fermi energy, the integrand contributes significantly only for — 1 < u < 1/VD. Therefore one can replace 
A2(u + uD) by its asymptotic value l/(— uD) m , and the Fermi function by its T — values up to terms of order 1/D. 
Y scales as 



D-^-Y cx D 1 '™. 



dD 

Only the skeleton diagrams of second order (m = 1) are relevant 



d f° 
D— Y otV 2 p'{u)du = V 2 p(0). 



Thus the second scaling equation is still valid and depends on the density of states only through its value at the Fermi 
energy. 

Magnetic impurity in a superconductor 

If the density of states of the host has still one energy scale D, but the jiensity at the Fermi surface vanishes 
according to a power law, this models magnetic impurities in superconductors^ with gap nodes: 



p(,) = (l + r)pg) r , (26) 



D 

p(oj)duj = 2Dp. 

D 



This modifies the scaling law (12) accordingly: Every skeleton diagram depends explicitly on D via the integration 
boundaries, which sorts out the NCA diagrams as the only relevant diagrams. In addition, every skeleton diagram 
has one factor D~ r per loop. Hence, the second scaling equation is 

D w Ff = {1+ r)(1 " l / N W-^ F f 

- r ~p y2 dW^ F f-^ + r yp v2 - ( 2? ) 
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The first scaling equation remains unchanged. In the limit of large N, this scaling equation is consistent with the 
results of FradkinEil, 

TK = D ( l -TTT7^)' r - ^ 



VIII. ALTERNATIVE DIAGRAMMATIC METHODS 

What conclusions can be drawn from this theory as to the reliability of other diagrammatic approaches, especially 
their ability to describe the crossover? 

A. NCA and the 1/N expansion 

If the system has N internal-degrees of freedom and exhibiiaJpcal Fermi-liquid behavior, the 1/N expansion becomes 
exact in the limit of large Nc3. Contrary to common belieffltJ, the 1/N expansion is not suited for the perturbative 
regime of high temperatures or magnetic fields because it fails to reproduce the diagrammatic renormalization-group 
results. That is so because the Kondo temperature itself is a function of 1/N. Hence it is not possible to describe the 
crossover within a finite-order 1 /N expansion. . .. . 

The NCA was considered so far as a "self-consistent" 1/N expansionOtJ. However, this does not explain why 
the NCA values for the static-magnetic susceptibility relative to their T = value agree so well with the respective 
renormalization-group resultsc3, even for N — 2. In view of Eq. (jl7|), this now becomes clear. 

B. Higher order skeleton diagrams 

One may ask whether it is possible to improve the NCA substantially by incorporating in T skeleton diagrams of 
higher order. Howeswer, one has to be aware that then the numerical problems in solving the self-consistency equations 
become formidablecJ. To date one has not succeeded to calculate x(0), and hence the Wilson number up to order 
1/N 2 by this diagrammatic approach. 

C. Non-self-consistent methods 

Obg way of incorporating higher-order skeleton diagrams in a theory for impurity systems was put forward by 
SasotJ in his T-matrix approach. The impurity propagators Rf were calculated within the NCA. For an observable, 
these NCA propagators were inserted into the respective skeleton diagrams of orders higher than 2. Such an approach 
cannot be derived from a variational principal. Hence it will not correctly describe the universal behavior of the 
impurity. 

D. Coleman's diagram technique 

In Coleman's approachEl to the Hamiltonian ([!]), the constraint rif + rib — 1 is enforced by adding A(n/ + rib) to H 
and sending A to infinity; consequently, Hf in Eq. (|l]) changes to 

H f = (e f + \)J2tiJ m + ^ + b . (29) 



The scaling equations ([l0|) and ( |l2| ) are rederived with the help of this diagram technique, and then Coleman's 
technique is used to discuss the conserving slave-boson approach of Kroha et a&B. 

Now that the unperturbed part H c + Hf is a one-particle Hamiltonian, the standard Matsubara-perturbation theory 
can be developed. The naked propagators from which the diagram technique is built are given by 
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R^(icj n ) = l/(ico n -e f -X), 
R§\iuj n ), = l/(iv p - A) 



-de. 



—D luJ n 



(30) 



Here ioj n — 2-n(n + 1/2)T and iv p = 2imT are the fermionic and bosonic Matsubara frequencies, respectively. The 
vertices are displayed in Fig. [I]. Again, a dashed line denotes Rm\ a wavy line -Rq°\ and a solid line represent the 
propagator of a conduction electron i?cm- This propagator already contains the sum over all momenta because of the 
local interaction with the impurity. 

Skeleton diagrams can be defined as above as containing no self-energy insertions. The variational principle follows 
with the help of the functional 



T 

T 



£ 



+ J^[ln {iujn - 
mn 

+ ]T{ln(k 2 



T,l^>(iuj n )R m (iu! n ) 



£/ — A - E m (iui n )) - ln(iw„ - e/ - A)] 



S cm (iw„)) - In 



R. 



(0) 



£ 

ji.r 



In {iu p — A — Eo(ifp)) — ln(zi/ p — A) 



(31) 



The factor 3 — 2/v arises because every term E^i? contains 3w/2 propagators. It is straightforward to show that T 
is stationary with respect to variations of the propagators iff the Dyson equations hold. By rearranging the internal 
summation-frequencies one has the following identity: 



E>m(iVn)Rm(iv n ) = £ cm (iw„)i? cm (ia;„) 

m,n m,n 

= (-1) y S (ii/p)-Ro(ii/p) 



(32) 



the minus sign in the second line corning from the additional fermion loopEl. 

As in Luttinger's original approachQ, one shows that T is equal to the difference of the interacting and noninteracting 
free energy F: 



T = F(V) - F(0) =: AF. 



(33) 



The Hilbert space is a sum of eigenspaces of the number operator n = J2 m fmfm + ■ The partition function 
can then be represented as a sum over these subspaces. The subspace n = describes the noninteracting system. In 
Coleman's original approach, the physical subspace n = 1 is projected out by 



d -AF 
urn — — 

a^oo de-' ix T 



lim Zf{X) = Zf. 

A — >oo 



(34) 



The limit is approached smoothly so that 



(35) 



1. First scaling equation 



Analogously to Sec. Ill A, the explicit derivatives with respect to pV , £/, and A are 
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pv r < 



= y^-Rm(^n), 



dpV 2 ~ 
d AF 



(36) 
(37) 
(38) 



To determine the T dependence of T, one has to bear in mind that each Matsubara frequency is proportional to T, 
that each summation over Matsubara frequencies gives a factor T, and that in a diagram of order v there are v/2 such 

summations. Furthermore, R^l depends on T and D only via D/T. Together with Eqs. ([3^), (|36|), (|37j), and (|3§| ) it 
follows that 



By using Eqs. ([54]) and (|35|), the analog of the scaling equation (|T(]) follows: 

d 



m 9 d T . 2 <9 



<9L> 



Z f = 0. 



(39) 



(40) 



2. Second scaling equation 



The functional T depends on the cutoff D only via the integration boundaries ±D in Rim- Hence 



D 



d AF 



E 

ran 



(-2) 



D 



D 



D iu> n — D 



x[R cm (iuj n ) - R^iun)]. 



The propagator difference R, 

Y 



Rem can be expanded as 

(-2) 



y pMy 

cm cm cm 



"™ "'cm I "cm 

In either case discussed below it will turn out that only the second-order term 

V 2 pT 



(41) 



(42) 



4mK 



N 



^RmiiUn + ilSp)Ro(ilSp) 



Ro,m(z) — 



is relevant in the universal limit of large D. The propagators i?o.m have the spectral decomposition 

dx . . 
rA),m(x) 

Z — X — A 

and are centeredll3 around A. After performing Matsubara summations, this gives 

d v 2 p 



(43) 



D 



dD T 



N 



-D0 £ 



dx dy pm(x)p (y) 

rn.a— ±1 1 

[/(x + A)+% + A)] [f[^-/M 
crD — x + y 



(44) 



with / and b the Fermi and Bose functions, respectively. The coefficient of e _/3A in Eg. ([44]) for large D, when 
combined with Eqs. (^) and ( p8|) yields because of Eq. ( pi] ) the desired scaling equation (|lg|). 

In Eq. ( p2| ) higher-order terms in £ are of high er order in e _,3A , and vanish relatively to S cm in the limit of large 
A. By the same arguments as given in Sec. IIIB, one can show that skeleton diagrams of higher order in E cm are 
irrelevant: Higher-order skeleton diagrams in Eq. ([42) contain at least three impurity propagators which have to be 
evaluated at iui n = ±D to make a contribution oc ijl) 2 , or they are of higher order in e _ ^\ 
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E. Conserving slave-boson approach 



Kroha et al.c3 developed, for the Anderson model, a method quite similar to the NCA. It is constructed from 
skeleton diagrams within Coleman's slave-boson technique. Kroha et ai, however, imposed the constraint only in 
the average (rif + n^n = 1. Their method then does not violate the gauge symmetries of the model like the usual 
slave-boson mean-field approachESl. The boundary condition for the free energy is 



d_ 



AF = (n f +n b ) H = 1. 



(45) 



This includes unphysical states, 
equation (p9[) now reads 



Consequently, the conduction electrons acquire a self-energy. The first scaling 



AF = 



m 3 <9 ir? 9 

T df + ^ f +pv dpv^ 



(46) 



In order to obtain any observable, AF has to be differentiated once more with respect to an external field. Hence this 
equation is the analog to Eq. (flO|). 

The second scaling equation can be derived from Eq. ( fii]) in the limit of large D. First it is summed over a. Only 
the term oc f(o~D) survive because the propagators R have a spectral weight centered around a finite value of A, and 
J2 a lim D ->oo D/(aD -x + y)=0. With Eqs. ©, (||), and @, one has 



D 



d AF< 2 ) -V 2 p/3 



dD T 



N 



E 



dx dy p m (x)p (y) 



x [f(x + A) + b(y + A)] 



(1 



Hn)v*p4-^-v>pP. 



de 



f 



T 



More generally, because X) CT nm -D^°o D/(iuj n + aD) = in every propagator of Eq. (^Tj) the Matsubara frequency 
iuj n can be replaced by aD with an overall factor f(o~D). In the magnetic limit one has A > |e/| to enforceEj both 
n/ w 1 and rib ~ 0. Therefore higher order terms in Eq. (|4|) are irrelevant, because they are either o(l/£>), or of 
higher order in e~^\ 

Hence the approximation also preserves the exact scaling laws. This explains why their results for the / propagator 
are so similar to the respective NCA resultsEJ. Moreover, one can now predict that taking into account more skeleton 
diagrams will not alter the picture qualitatively. 



IX. SUMMARY 



The Anderson-impurity model shows, in a nutshell, the difficulties when dealing with strong correlations. One 
encounters the same problems as in high-energy physics: The perturbation theory of observables at energies T in the 
coupling constant diverges logarithmically, oc In D/T. The limit of zero energy T in solid-state physics corresponds to 
the limit of infinite cutoff D in high-energy physics. If the Hamiltonian can be shown to be renormalizable, it means 
that there are finitely many energy scales in the system, at least in the perturbative region. If there is but one scale, 
the system can be described by a "running" coupling constant J(D), expressed with the help of the (3 function: 

Reducing the cutoff D of the system does not change its physics as long as the effective coupling J(D) is changed 
according to Eq. (^7|), to keep the energy scale constant. However, outside the energy regime where the perturbation 
theory is valid it has never been shown that such a [3 function exists at all. 

In this paper, the variational functional of Luttinger and Ward was used to prove both the existence and form of 
the [3 function, or, in other words, the exact energy scales of an impurity-system were determined. It turned out that 
neglecting the vertex corrections (NCA) suffices to proof universal behavior, and that families of higher-order skeleton 
diagrams do not alter the energy scales. 
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The method differs therefore from the conventional renormalization-group approach. Within that approach, the 
fixed points of a flow of effective Hamiltonians is studied perturbatively as the energy scales of the system are varied. 
The crossover in inaccessible by that method. 

In comparison to that, the skeleton diagrams can describe the crossover very well, as was demonstrated. However, 
the theory cannot predict the nature of the fixed-point Hamiltonian. 

Further work on this subject will concentrate on the question of whether the NCA can also be justified as a means 
to solve the effective impurity model onto which the infinite-dimensional versions of correlated electron systems can 
be mappeoOCj. In particular one would like to learn whether the periodic Anderson Hamiltonian exhibits heavy 
fermion behavior in this limit, or if this is only an artifact of the approximations used. Also, it should be within the 
reach of the theory to decide which class of diagrams should be used to describe the problem of two impurities in a 
metal. 
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APPENDIX A: IRRELEVANCE OF HIGHER-ORDER SKELETON DIAGRAMS 

The proof uses the spectral decomposition of the ionic propagators Rf. Every skeleton diagram depends explicitly 
on D only via the integration boundaries ±D coming from integration over the conduction-electron lines, as shown 
in Eq. (ph. For a = 0, . . . , N, let ^ R a be such a skeleton diagram of 2nth order with p closed fermion loops and 
n ionic propagators Rf. 

dZ -e-^\z)R a {z) 



2m 

p ny2n 

N n -P T 2m 



£-e~" z n / /fe)^- ipm * + 5> ■ (ai) 

J=l D *=1 \ j } 



Here J dej denotes the integration over the jth conduction-electron line. Furthermore, the term m the argument 
of Ri indicates that the sum runs over j if the ith propagator sits under the jth conduction-electron line, as in 
Fig. ||, the numbers indicating the €j. The ith propagator Ri of the respective skeleton diagram has the spectral 
decomposition 



(») roc / / 

Ri{z + J2^) = J d\iPi(\i) j I 



For that reason, the line integral gives only contributions from the real poles 

(fc) 

z = uj k = A fe - J^ej. 



For real u, there holds the relation 

/.oo W-X 

where the Cauchy-principal value of the integral has to be taken. Hence the fcth pole of the line integrals gives 



™ r D r°° a , , 

l[ / f{e 3 )de 3 / du e - p{uJ -^ €i) p k (uj) 

■_! J —D J — oo 



J'=l 

(fc) (i) 
i^k j j 
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Because of /(ej)e^ £j = /(— e,-), one shifts those variables ej to —ej, which occur in but not in 

where now runs over j if the propagator Ri sits under the jth conduction-electron line in the respective skeleton 
diagram with cyclic permutated vertices. In such a diagram, the propagator Rk is the outer one. 

Exactly here one uses the fact that only whole families of skeleton diagrams are considered: It was just shown that 
for every family F of 2nth order skeleton diagrams their contribution to the line integral in Eq. (||) is 



Tr f j ^ e -0^ n ' F \z)R f {z) = 2njdw 



Here, the operator Re is defined to replace every propagator in its argument by its real part. The factor 2n arises 
because every skeleton diagram occurs in a family 2n times. In T, every skeleton diagram of order 2n depends 
explicitly on D only via the integration boundaries ±D of its n integrations over the conduction-electron lines. The 
ith integral over a conduction-electron line can be written as 




D 
-D 



Re£< 2n ' F) M= I / I Ai(w) / deJ(e l )A 2 (iu + e l ), (A2) 



where the part Ai of the respective diagram does not depend on e^, that is, its ionic propagators are lying on the left 
or the right of the ith conduction-electron line. The propagators of part A2 lie under the ith conduction-electron line. 
(/)" hints at the other n— 1 integrations, weighted with the respective Fermi functions. One differentiates with 
respect to D by evaluating the integral over the ith conduction-electron line at its integration boundaries ±_D, 



2n I due~ l3u p 0>m (cj)[A 1 {uj)Df{aD)A 2 (u} + aD)] 

~ LI •* 



CT=±1 

The integral is to be weighted with 1/Zf in order to yield the D derivative of the respective contribution to the free 
energy. Because e^ E °Zf goes for low temperatures to 1, the spectral densities po,m are weighted with e^^ E °~ u \ For 
low temperatures T <C D, this weighted spectral density contributes to the integral only for frequencies u) < Eq. 
Because the other integrations over conduction-electron lines Cj are weighted with Fermi functions f{ej), respectively, 
for low temperatures T <§C D the real part of a propagator $tRf(u> — D + J2j e j) contributes in A 2 (w — D) only for 
frequencies 

oj-D<E -D< -D. 

In this frequency interval one can replace by its naked counterpart l/(u> — Hf — D), and hence can estimate 
its contribution to Ff from above as 1/D, The term oc f(D) does not contribute for T -C D. It therefore can be 
estimated as 

LD±l¥»R f ocD/iy«. 

Hence at least one must equal 1, otherwise the skeleton diagram is irrelevant. However, such a diagram is a 
skeleton only if it is of second order; otherwise it would have a self-energy insertion. This was shown. Eq. (113) follows. 



APPENDIX B: NCA AT ZERO TEMPERATURE 



1. NCA differential equations 



The NCA equations (||) together-svith the Dyson equations ([I|) constitute a self-consistent system of equations 
which have been solved numericallyE3. For T = the NCA integral equations (||) in a magnetic field are transformed 
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into differential equations by substituting for the Fermi function a step functionOI§li-S, by introducing the negative, 
inverse of the propagators i?o, m , 



Y (z) = ^\z)-z, 
Y m (z) = Y, ( ${z)+e f -z, 

and Y m = Y m + rngpgh for —j < m < j = (N — l)/2, 

d Yi _ 1 P y2 y _ 1 

da; ° N m Y m + mgp B h 

d - 

The inverse propagators have the asymptotic forms 



(Bl) 



N 

pV 2 
N 



Y; 



Y m (uj) « -u> + e f 
Y (uj) M -U! 

for large, negative u>. The NCA differential equations have, up to terms of order 0(1/ D) the integral 



(B2) 



(B3) 



Y 
pV 2 



Y n 



Y, 



pV 2 N 



P V 2 

*2> 



Y m + rngpsh 



f = *t+[ 1 -jf)P V ln 



P V 2 

D 
~pV 2 ' 



pV 2 



(B4) 



The value of follows by inserting Eq. (fB3|). Terms of order 0(1/ D) are neglected because the universal behavior 
of the system is investigated. The integral of the NCA differential equations contains the energy scales pV 2 and Tr 
viaEd ln(T^ / pV 2 ) — ej/pV 2 . It is a nontrivial task to solve this differential equations numerically in the universal, 
magnetic limit 

D — > OO, 6/ — > — oo 

e f V 2 
lim = 0, i — : = J = const.. 
Z^oo D | e/ | 



2. What is the ground-state energy of the NCA? 



The ground-state energy of the NCA is defined as 



£ NCA = lim Ff CA . 



(B5) 



However, one has not yet succeeded in deriving an-expression for Eq CA via that route, but merely solves the NCA- 
differential equations (BS). It has been conjectureda that the lowest, real zeros of the inverse propagators lo,m define 
the NCA ground-state energy. This is now proved: Because of Eq. (|B5|) Eq CA fulfills 



1 = lim e pE ° 

T->0 



-0ui 



p? ca h+E^ ca (^ 



duj. 



Hence, for T — there exist the (positive) spectral densities 

P0,m M = e ° P0,m M, 



(B6) 
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and therefore the spectral densities Pq^ a van i sn f° r T — and ui < Eq CA . This means that the inverse propagators 
io,m are real for those frequencies. In addition, p^^ vanish for w > Eq, because Pq ^ A remain finite. However, from 
the existence theorem for solutions of the differential equations (B2) it follows that there exists at least one value ujq 
such that the following hold. 

(1) lo,m h ave a zero in loq and become complex above, if they are real for large, negative u> — which they are in 
view of Eq. (|B^). Because it was just shown that the Yb.m are real for u> < Eq CA , they are real for u> < loq as well. 

(2) vanish for uj < u) , therefore w < Eq GA . 

(3) Po,m A vanish for uj < uj , therefore ujq > Eq CA . 

To sum up: ujq is the lowest common zero for Yq m , and at the same time the NCA ground-state energy. 



3. Parametrization of the NCA ground-state energy 



Hence the well-known expression for the common zero of F .m^""'" can be used as the NCA ground-state energy 
and is parametrized in the following manner: Define the W- function for positive x as 



W(x) exp = x for x > 0, 



with asymptotic behaviors 



W(x) = x + o(x 2 ) for 



+ , 



W'[cxp(.T;)] = x — ln(.x') + o(l) for x — > 30. 



The integral (B4) can be solved for Yq because both Yq and Y m are positive for ui < Eq CA : 



NV- 2 Y Q 
= W ( Nexp 



rv *\ N , Y^i Ym + mg/jBh 



pV 2 



(B7) 



(B8) 



(B9) 



can now be parametrized with the help of 



E = D 



dw = D 



D 



jgfj, B h 



Y m (-D) 



dw 



-dY r 



the NCA differential equations (B2), and Eq. (JB9|) as 



E Q = ef- jgpBh - 



^5^ 



dw 



1 + W \ Nexp 



pV 2 



(BIO) 



For vanishing magnetic field there exists another parametrization because Y m can be written as a function of Yq with 
the help of the W function: 

r-D 



Eq = 



did 



1 + W [exp 



I L *\ 1 1 1 ^ 

[uj + tf)—rp; + — In ■ 



(Bll) 



f, pV 2 ' N^pV 2 _ 

The NCA ground-state energy fulfills the equations ( |Io| ) and (^2|), as can be checked with the help of (Bll). The 
ground-state energy is not an universal function because of the constant term — pV 2 in Eq. (|l2|). The ground-state 
energy up to order 1/N and 1/D follows asa: 

Eq = e f + pV 2 W -prcxp 



+ 



D 
~p~V 2 

W(y) 



N J [1 + W{y)Y D 
D 

W 2 ^ V pv 2 



In dx + 0(1/ D) + 0{1/N 2 ) 



(B12) 
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4. An analytical expression for the static magnetic susceptibility of the NCA 



a. Universality 



With the help of the Kondo temperature (|l4|), Eq. ( |B10 ) reads, in the magnetic limit Tk / pV 2 — > after the 
substitution lu /Tk — x, 



E = e f - jgpBh -T K 

J J: 



(D+e f )/T K 

gfi B h/T K 
1 



dx 



(B13) 



1 + W { N exp 



Em. 



After differentiating with respect to the magnetic field h the integrand decays rapidly enough so that one can replace 
the upper limit of integration by oo. Hence the magnetization scales as a function of h as 



M(h,D,e f ,V,N) = M 



(B14) 



This shows explicitly that already the skeleton diagrams of second order give the exact energy scale for the spin 
degrees of freedom. 



b. Small magnetic fields 



Using the parametrization (BIO), the magnetization M(h) — —dhEg vanishes for h = because of J2 m m ~ 0. ^ ne 
second derivative gives the static magnetic susceptibility for h — 0. With the abbreviationMl 



\n]N = ( 3 /^) 2 ^m 2 , 



from Eq. ( B10| ) one has, after substituting NYo/pV 2 = W , using the NCA differential equations, and finally partially 
integrating, 

x(o)=v— r 2w(y)+i dx 

XW 3%Wo W(y)[l + W(y)] 3 

y=(^^p(^)- (bis) 



\pV 2 J pV 2 \pV 

In the magnetic limit, y is very small for x < \e*f\ because of — ^> pV 2 , and above that W > 1 and the integrand is 
small. Therefore one can replace the W function by its argument, and extend the integration to oo: 



X NCA 



= ^u 2 ±r(l-l/N). (B16) 
A 1 K 

The exact result 

^ = l^WTW) (B17) 

was obtained by fitting the results of the Bcthc ansatz to perturbation theory in the nonmagnetic limiti. Both results 
coincide up to order lJJf, because the NCA contains all diagrams up to that order. This contradicts the claim of 
Kuramoto and KojimaEj that the NCA would yield the exact result for x(0) in the magnetic limit. 
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APPENDIX C: COQBLIN-SCHRIEFFER MODEL 



1. Variational functional 

The Schrieffer- Wolff transformatior0 projects - up to a constant - the Anderson model onto the Coqblin-Schrieffer 
model in the magnetic limit e/ — > — oo and constant J = V 2 /\ef\, 



H^H CS ~2 



pJD 

N 



(CI) 



]J r „ = V e c+ c +— Vr+r f+f 

/ t ^P^pm^pm 1 / { ^pm^qnJn J n 



M<-D,' 



pq,mn 



For that the spectrum of H is-shifted at (— e/) by shifting the argument z in T at e/, and the propagators and 
self-energies are transformed astJ 

— Z (z + e f ) ^U{z), 
e f 



1 



Ro(z) 



l + f--±Z ( z + e f ) i-n(z)' 

R m (z + £/)—> R m {z) = ^ , r . 

A 2nth order diagram carries the prefactor (— l) n (pj) n . Because |e/| ^> D,\z\, the term z/ej was neglected in 
CfRo(z + €f) and such the charge degrees of freedom of the impurity are projected out. The variational functional T 
now reads 



T 



-»/s-*{E(i-i)=& , w«.c 



-In 



-^S^(z) +ln l-£lT»(z) 



(C2) 



The saddle-point property of T is shown as Bickers byEf 



2. Skeleton diagrams of second order 



The analog to the NCA is called the "self-coasistent ladder approximation" 0. The NCA equations follow from 



Eq. (H) after projecting as in the last paragrapl 



raprM 



pj 



E m (z) = -^/ /(e 



D 



1 



D i-n(z + e) 

D 



dc 



(C3) 



n W = -77E J D f(e)R m (z + e)de. 



The first scaling equation follows as in Sec. [II A 



d d 
F f = T——Ff + D— — Ft. 



The asymptotic behavior of the self-energies is estimated from Eq. ( |C3| ) for \u>\ -C D as 

S m (w - D) = -pJD/N + 0{J 2 ) 
H(u-D) = 0(J InD). 



(C4) 



(C5) 
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Because of Eq. (CI), z has to be shifted in every propagator by —2pJD/N, to describe the Coqblin-Schrieffer model. 
Therefore 

R m {uj-D)= ^—i— ^ forH«A (C6) 
Ro(w - D) = 1 for|w| < L>. 
Higher orders in J are irrelevant as will be shown below. It follows within the NCA that 

D *Z, = -PJ I —e-^ ( Ro( - Z ] DR m (z) 
dD S H J 2m \l + pJ/N my 1 

Hence the scaling equation for the impurity part of the free energy is 

n ^ fNCA _ {PJ) 2 ® F NCA PJ n r r7 \ 

°dD F f ~ 1 + pJ/Ndpl f ( } 



This equation does not change in the universal limit J D if higher orders of J in Eq. (C6) are taken into account. 
The following scaling law holds therefore for the impurity part of the free energy: 

(F f ~E f CA (T,D,J) = (F f -E f CA (J^j (C8) 

with the Kondo temperature of the Coqblin-Schrieffer modeli 

yes = D \fp~j cxp [—1/pJ]. (C9) 

For other observables, one has to couple H to suitable external fields. In particular, the Kondo resonance is reproduced 
qualitatively correctly. 



3. Skeleton diagrams of higher order 



The two limits D — > oo an d e t — > — oo are not interchangeable, as can be seen from comparing the Kondo 
temperatures of Eqs. (|lj) and (C9). Skeleton diagrams of higher order than two are relevant in T. The reason for 
that is the asymptotic behavior of Rq which goes to one at the cutoff. If the diagram in Fig. || is logarithmically 
differentiated with respect to D, the contribution of the second co nduc tion-electron line does not vanish for large D. 
In particular it is not possible to prove now that the energy scale (C9) is the exact one by considering only skeleton 
diagrams of second order. In fact, the NCA still predicts for N = 1 a low-energy scale (C9) although there is none. 



4. NCA at zero temperature 



The derivations are as in Appendix |[ The inverse of the pseudopropagator n is defined to be Yq = 1 — n. The 
NCA differential equations are up to terms oc 1/-D, 



9 Y _ ~PJ 



(CIO) 



d_ 



Y 

J IT, 



PJ, 



-{l-pJ/N)- H —Y-\u). 



With Y m = Y m + mgpsh and D = D(l — pJD/N) the NCA differential equations have up to terms of order X/D the 
integral 



1 1 



[—--)+ -W, 



pj N 



1 



N 



N ^ 



Y m + mgp, B h 



D 



1 1 
~p~J~N 



(Cll) 
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The ground-state energy is expressible as 



E = -A+ I dw, 

A 



where A is a still arbitrary constant. If J -c A <C D and \\mr>^ca A = oo, the integral (Cll) can be used as in 
Appendix |b] to yield, in the universal limit of small p J, 



,NCA^_ 1..2 i r(i-i/jv) 



* ( ° ) = i^# V~e • (C12) 
which is up to 0(1/N 2 ) identical! to the result of Rasul and Hewson, where T(l — 1/N) is replaced by l/r(l + 1/N). 
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FIG. 1. Vertices for the Anderson model. A dashed line represents the naked propagator of the occupied f 1 configuration 
with internal quantum number m. A wavy line represents the naked propagator of the unoccupied /° configuration. A solid line 
represents the propagator of a conduction electron with internal quantum number m. Every diagram has a spine of alternating 
wavy and dashed impurity propagators. The conduction-electron propagators carry no self-energytJ. 

FIG. 2. Skeleton diagram for E^, 2) (R m (z))R (z). The last vertex can be identified with the first because of the trace over 
the / configurations in Eq. Q. The double, dashed and double, wavy lines represent dressed propagators. 
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FIG. 3. Skeleton diagram of order 6 for Eq-Ro 
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